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Abstract. According to a currently popular paradigm, 
nuclear activity in quasars is sustained via accretion 
of material onto super-massive black holes located at 
the quasar nuclei. A useful tracer of the gravitational 
field in the vicinity of such central black holes is avail- 
able in the form of extremely dense gas clouds within 
the broad emission-line region (BLR) on the scale of 
~ 1 parsec. Likewise, the radio sizes of the lobe-dominated 
radio sources are believed to provide a useful statisti- 
cal indicator of their ages. Using two homogeneously ob- 
served (and processed) sets of lobe-dominated radio-loud 
quasars, taken from literature, we show that a positive cor- 
relation exists between the radio sizes of the quasars and 
the widths of their broad Hp emission lines, and this cor- 
relation is found to be significantly stronger than the other 
well known correlations involving radio size. This statisti- 
cal correlation is shown to be consistent with the largest 
(and, hence, very possibly the oldest) radio sources har- 
boring typically an order-of-magnitude more massive cen- 
tral engines, as compared to the physically smaller and, 
hence, probably much younger radio sources. This infer- 
ence is basically in accord with the " accreting central en- 
gine" picture for the radio-loud quasars. 

Key words: Line: profiles- Galaxies: active-Galaxies: 
nuclei -quasars: emission lines-quasars: general-radio con- 
tinuum: galaxies 



1. Introduction 

Being too compact to be resolved with even the most ad- 
vanced optical telescopes, the structure and kinematics of 
the broad emission line region (BLR), a prime feature of 
quasars, continues to be a major enigma in the AGN re- 
search (see, e.g., Brotherton, 1996; Marziani et al., 1996; 



Corbin, 1992). Nonetheless, it is a key ingredient to the 
theoretical models that seek to explain various observable 
properties of quasars, e.g., their intense 7-ray emission 
(e.g., Dermer & Schlickeiser, 1993; Ghisellini & Madau, 
1996). Deciphering the BLR geometry has therefore been 
a key objective of many observational programmes. One 
strategy for this is the so called 'reverberation mapping' 
(e.g., Peterson, 1993; Koratkar & Gaskell, 1991), from 
which a positive dependence of the BLR size on the bolo- 
metric luminosity has been inferred: r ~ 0.06L4g 5 pc (Net- 
zer & Laor 1993), where is the luminosity expressed 
in the units of 10 i6 erg.s . 
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An early indication about the BLR geometry came 
from an empirical study of a heterogeneously selected sam- 
ple of radio-loud quasars with the characteristic core-lobe 
type radio structure. The study revealed a statistically 
significant anti-correlation between the prominence of the 
radio core relative to the lobe, and the FWHM of the 
HP broad emission line (W). It was thus inferred that 
the BLR clouds are predominantly confined to a rotating 
disk-shaped region surrounding the quasar nucleus and 
oriented roughly perpendicular to the jet axis (Wills & 
Browne, 1986; hereafter WB86). The significance of this 
correlation was found to be considerably lower in a sub- 
sequent study, however (Jackson & Browne, 1991b). In a 
recent work, it has been proposed that the absolute vi- 
sual magnitude, M v , of the quasar (plus its host galaxy) 
provides a more reliable measure of the intrinsic power of 
the central engine, and therefore the beamed radio core 
flux normalized by M v is a better indicator of the ori- 
entation of the jet relative to the line-of-sight (Wills & 
Brotherton, 1995; Brotherton, 1996). Adopting this new 
parameter for the core-prominence and employing a larger 
sample of quasars , these authors have found a conspicu- 
ous anti-correlation between the core-prominence and the 
W(HP), thereby supporting the disk-like BLR geometry 
inferred earlier by WB86 (at least for the Hp emitting 
clouds). 
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Table 1: Sample of lobe-dominated quasars (42 quasars) 
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It is now widely believed that the large widths of the 
BLR emission lines are a manifestation of the deep grav- 
itational potential at the centers of quasars, possibly due 
to super-massive black-holes. The accretion of the mate- 
rial postulated for sustaining the quasar luminosity is ex- 
pected to steadily increase the mass of the central engine 
during the lifespan of the nuclear activity. If this conjec- 



ture is basically right, the question arises: do we see in the 
data any evidence for the postulated growth of the cen- 
tral mass (e.g, in the dynamics of the BLR clouds) ? In 
the present study we shall examine this issue by employ- 
ing published measurements of radio-loud, steep-spectrum 
quasars. While the needed reliable estimator of the age is 
generally not available for individual sources, the overall 
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radio size can serve as a useful statistical measure of age 
(since most radio sources are believed to steadily grow in 
size with time; cf. Sect. 3). Thus, we wish to examine here 
the nature of any relationship between the observed lin- 
ear sizes of steep-spectrum radio quasars and the widths 
of their broad H/3 emission lines. 

2. The datasets of radio-loud quasars 

In order to minimize the projection effects on the mea- 
sured radio sizes, we confine our study to lobe-dominated 
quasars (LDQs). Measurement of the other parameter, 
namely, W, the FWHM of the broad H/3 emission line, is 
rather complicated, due to the contamination arising from 
the lines of Fe II and the narrow component of H/3. There- 
fore, adopting uniform observing strategy and profile ex- 
traction procedure for the entire sample is an important 
pre-requisite for a meaningful interpretation of the data. 
Using two such datasets on H/3 line widths available in the 
literature (Brotherton, 1996; Jackson & Browne, 1991), for 
which we could also find the requisite radio data, we inves- 
tigate the relationship between different radio properties 
and W(H/3) (Table 1). The largest linear sizes, I, of ra- 
dio emission associated with these quasars are taken from 
published radio maps (Table 1); these correspond to the 
cosmological parameters: H a = 50 km s _1 Mpc -1 and 
q = 0, as adopted throughout this paper. Details of the 
two sets of LDQs (all having M v < —23 ) are: 

(a) B(96): This first set is derived from the sample of 
Brotherton (1996), which itself was based on the QSO 
compilation by Veron-Cetty & Veron (1988). His selection 
criteria were: (i) core-fraction, f c should be available in 
the literature (f c is defined as the ratio of the core to 
extended flux density, measured at 5 GHz in the rest-frame 
of the quasar), (ii) V < 18 mag, (iii) z < 0.95 and (iv) 
declination > —20°. For the selected 60 quasars, optical 
spectra were taken with a typical spectral resolution of 
2.5 A. Out of these, we have selected all the 31 LDQs 
(log/ c < 0) having an absolute magnitude M v brighter 
than -23. 

(b) JB(91): This set is derived from Jackson & Browne 
(1991a, b), whose sample consists of low-resolution, (20 — 
25A) spectra of 53 radio-loud quasars selected using the 
following criteria: (i) log(R) > 1 [R is defined as the ratio 
of the flux densities at 6 cm and 4400 A, in the rest-frame], 
(ii) V <18-mag, (iii) z < 0.86, (iv) right ascension between 
16 and 13 hr, and (v) declination > —30°. The sample 
contains 23 LDQs with available H (3 profile measurements 
by Jackson & Browne (1991b) who determined W using 
two different profile extraction procedures. We adopt here 
the values of W obtained using the four Gaussian fitting 
method as it corrects for the the narrow line component. 

As seen from the last column of Table 1 , twelve quasars 
are common to the two datasets. The ratio of H/3 line- 
width determined by B(96) to that by JB(91), designated 
as £, has a mean value of 0.68±0.16. Thus, in most of 



the cases, £ is close to the mean value, the only exception 
being the quasar 0710+118 (see, Sect. 3). We form a com- 
bined dataset by merging the quasars present exclusively 
in the set JB(91) with the set B(96), after multiplying 
their quoted H/3 widths by the mean value of £ found 
above. 

Table 1 gives a consolidated list of all the LDQs used 
in our study. The information provided includes: the IAU 
name of the quasar, its redshift (taken from Veron-Cetty 
& Veron, 1993), the largest angular size (LAS) measured 
from the radio maps (references are cited in the footnote). 
Column (4) tabulates the value of log f c . Column (5) 
gives log R v , defined by Wills & Brotherton (1995) as 
logi?„ = log(L core /L opt ) = log(L core + M v /2.5) - 13.69, 
where the absolute visual magnitude, M v , is taken from 
Veron-Cetty & Veron (1993) and listed in column (6). 
Columns (7)&(8) tabulate the published widths of the 
broad H/3 emission line, taken from the two datasets dis- 
cussed above. Column (9) lists the value of £, the ratio of 
W{H(3) measured by B(96) to that by JB(91) for the 12 
quasars which are common to the two data sets. 

3. Results 

Fig. 1 shows the radio size, I, plotted against the parame- 
ters, Rv, f c , W(H/3) and z for our sample of 42 LDQs. The 
filled circles are the data from B(96) and the open circles 
refer to the data from JB(91). The values of W shown by 
the open circles have been rescaled, as discussed above. 
Table 2 shows the results of the non-parametric Spear- 
man rank correlation tests between the different parame- 
ters, for B(96) as well as the combined dataset. For each 
case, the upper line gives the correlation coefficients and 
the lower line gives the two-sided significance level of its 
deviation from zero (smaller value implies a stronger cor- 
relation). The remarkable trend noticed from Fig. 1 and 
Table 2 is that I is found to correlate more strongly with W 
than with any of the other parameters, namely f c , R v , and 
even redshift, z. The I — W(H(3) correlation thus emerges 
to be of the primary statistical significance , in agreement 
with the preliminary results reported by Gopal-Krishna 
(1995) (see, also, Gopal-Krishna & Srianand, 1998). It 
may be recalled that using a much more limited dataset, 
Miley & Miller (1979) had earlier noticed a positive cor- 
relation between I and W(H(3). However, almost <~ 40% 
of their sample was comprised of core-dominated quasars 
and hence, the linear sizes used by them are likely to be 
influenced by projection effects to a much greater degree 
(and, consequently, be rendered less suitable as a measure 
of source age, vis-a-vis the present study where only lobe- 
dominated quasars have been considered). Note that the 
quasar 0710+118 (marked with an arrow in Fig lc) shows 
a large deviation from the general trend. Since this quasar 
is the only case where W(H/3) measured by B(96) is larger 
than that found by JB(91) (Table 1), we have also plotted 
in Fig lc the scaled JB(91) value of its W(H/3) [the open 
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Fig. 1. (a-d) The plots of R v , fc, W(H/3), and z, versus the radio size (I) for our combined dataset of lobe-dominated quasars. 
The filled circle are the data from B(96) and the open circles represent the data from JB(91) (Sect. 2 and 3). 



Table 2: Spearman rank order correlation coefficients 



Sample 


Size 


l-z 


I - logifc) 


I — log(R v ) 


I - W 


log(fc) - W 


log(R v ) - W 


B(96) 


31 


-0.163 


-0.207 


-0.376 


0.771 


-0.254 


-0.397 






0.373 


0.254 


0.030 


2.4 x 10~ 7 


0.161 


0.024 


combined 


42 


-0.237 


-0.111 


-0.420 


0.765 


-0.183 


-0.324 


dataset 




0.127 


0.447 


0.005 


2.3 x 10" 9 


0.240 


0.034 



circle connected by a dotted line] and this point is fully 
consistent with the general trend. 

In order to understand possible implication of the I— W 
correlation, we first recall the general consensus that a vast 
majority of powerful extragalactic radio sources grow in 
linear size with increasing age (e.g., Rcadhead, 1995; Fanti 
et al., 1995). Thus, to a first order, the observed linear 



size, I (i.e., the separation between the outermost peaks 
in the radio lobe-pair) can be regarded as a meaningful 
statistical indicator of the age (e.g., Best et al., 1996). For 
individual sources, radio spectral gradient is also used as a 
measure of the age (sect. 4). However at present such data 
are largely confined to very extended radio sources with 
prominent diffuse lobes, mainly radio galaxies where the 
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contamination of the lobe emission by the jet component is 
usually quite small (e.g., Alexander & Leahy, 1984; Carilli 
et al., 1991; Liu & Pooley, 1992; Rawlings & Saunders, 

1991) . 

A potentially complicating factor in the use of I is the 
foreshortening caused due to the (unknown) projection ef- 
fects. However our selection criterion, which excludes core- 
dominated quasars, would effectively minimize the role of 
projection. It may be recalled that I is known to anti- 
correlate with z for quasars, as reported, e.g., by Wardle 
& Miley (1974), which is usually interpreted in terms of an 
increase in the density of the ambient medium. However, 
this may not play a major role in our data as it covers a 
much smaller range in redshift (Table 1). This expectation 
is supported by the apparent weakness of the I — z correla- 
tion in our data (Table 2). Radio size, I, is also known to be 
anti-corrclatcd with the radio core prominence, / c , which 
is consistent with the idea of beamed core radiation (e.g., 
Kapahi & Saikia, 1982; Browne & Perley, 1986; Lister et 
al., 1994). Again, this trend is at best only weakly present 
in our datasets (Fig. 1; Table 2), which is in accord with 
the expectation that projection effects should be of minor 
importance in the case of lobe-dominated quasars. This 
is further supported by the apparent weakness of I — R v 
anti-correlation in our datasets (Fig. 1; Table 2). 

4. Discussion 

An important inference from Fig. 1 is that over the lifetime 
of a typical powerful double radio source, the FWHM of 
the HP emission line from the BLR undergoes an increase 
by roughly a factor of 3-4. This broad line is presumed to 
arise from the vicinity of the central super-massive black- 
hole (SMBH) and its width is commonly attributed to the 
gravitational influence of the SMBH, in accordance with 
the standard paradigm for AGN (however, for radio-quiet 
AGN, alternative schemes have been put forward by, e.g., 
Tcrlcvich et al., 1992). In the past, widths of the BLR 
emission lines have often been used as the key tracer of 
the mass of the central engine of AGN (e.g., Dibai, 1984; 
Wandel & Yahil, 1985; Padovani & Rafanelli, 1988, Perry, 

1992) . 

Here it may be recalled that the radio spectral map- 
ping programmes have yielded typical lifetime of order 
10 7 yr for powerful double radio sources (e.g., Alexander 
& Leahy, 1987; Carilli et al., 1991; Liu & Pooley, 1992). 
However, independent considerations based on the dynam- 
ics of the radio hot-spots suggest that the true ages are 
probably an ordcr-of- magnitude larger (i.e., ~ 10 8 yr). 
This is borne out, for instance, from a recent analysis of 
the lobe-length asymmetry, in which the near-side radio 
lobe in individual sources could actually be identified us- 
ing the observed one-sided jet (Scheuer, 1995). Indeed, 
several authors have argued that the analyses of spectral 
gradients in radio lobes could substantially underestimate 
the ages, since they usually incorporate the simplistic as- 



sumption of a uniform magnetic field. As indicated by 
the well-resolved maps of radio lobes and jets, the field 
is likely to be concentrated in filaments occupying just a 
tiny fraction of the total volume of the lobes (e.g., Per- 
ley, Dreher & Cowan, 1984; Owen, Hardee & Cornwcll, 
1989). Consequently, the energetic particles would spend 
most of the time in the weak-field region, thereby reducing 
the synchrotron losses by a large factor (see, e.g., Eilek, 
Melrose & Walker, 1997 and references therein; Gopal- 
Krishna, 1980). Likewise, the spectral ages may also get 
underestimated if (large-scale) spatial gradients of mag- 
netic field are present within the radio lobes /jets (Wiita 
& Gopal-krishna, 1990). Based on all these considerations, 
it seems more likely that the true ages of powerful radio 
sources are close to 10 8 yr. For the case of Eddington ac- 
cretion rate one would then expect typically an order-of- 
magnitude increase in the mass of the SMBH during the 
lifetime of the radio sources, since the e-folding time for 
the black-hole mass is: 

t E dd = \MEddfM\~ 1 = (4.5 x 10 7 yr) (ry/0.1), (1) 

where ry is the radiative efficiency. We argue below that 
even though the persistence of an Eddington accretion rate 
is not supported by the data, the inference about an order- 
of-magnitude increase in the mass of the central engine 
(i.e. the region within the BLR) during the radio source 
lifetime probably still holds. 

In the event of Eddington accretion rate the increase 
in the mass of the central black-hole with time would re- 
sult in a similar increase in the luminosity. As a result, 
one would expect a positive correlation between L opt and 
Z, which is not evident in our combined data set (Fig. 
2a). This weakens considerably the case for a persistent 
accretion at the Eddington rate. On the other hand, in 
order to estimate the increase in the mass of the central 
engine one needs to take into account the result from the 
'reverberation mapping' of nearby AGNs, which has re- 
vealed a positive dependence of the BLR radius (r) on 
the luminosity (L): r ~ 0.06L^ 5 pc (Nctzer & Laor 1993), 
where is the bolomctric luminosity expressed in the 
units of 10 46 erg.s _1 . Combining this with the relation 
v 2 oc GM/r gives: M oc v 2 ^L 46 . Fig. 2b shows a plot 
of M versus linear size, I which is a measure of age of 
the quasar (note that L has been approximated by opti- 
cal luminosity, cf. Joly et al., 1985). A positive correlation 
between M and I at a 4.5<r level is seen in our combined 
dataset of LDQs. This suggests that over the lifetime of 
a typical radio quasar, the mass of the region inside the 
BLR does go up by about an order-of-magnitudc. These 
estimates could be further sharpened when the estimates 
of bolometric luminosity for the individual quasars be- 
come available. Another interesting improvement to this 
study would be to induct additional quasars with very 
large radio structures. That would shed light on the late 
evolutionary stages of the central engine of quasars. 
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Fig. 2. (a,b): Plots of radio size (I) versus the optical luminosity and the mass of the central engine, respectively, for our 
combined dataset (see Sect. 4). The filled circle are the data from B(96) while the open circles refer to the data from JB(91), 
as explained in Sect. 2. 
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